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1 .  INTRODUCTION 


The  problem  of  testing  the  hypothesis  of  the  equality  of  the  mean 
vectors  of  several  multivariate  populations  with  a  common  covariance 
matrix  received  considerable  attention  in  the  literature.  The  test 
procedures  are  based  upon  certain  functions  of  the  eigenvalues  of  the 
multivariate  analysis  of  variance  (MANOVA)  matrix.  In  the  univariate 
case,  the  MANOVA  matrix  reduces  to  the  ratio  of  the  between  group  and 
within  group  sums  of  squares.  The  joint  distribution  of  the  eigenvalues 
of  the  MANOVA  matrix  in  the  noncentral  case  is  useful  in  studying  the 
power  of  the  tests  for  the  equality  of  the  mean  vectors.  This  distribu¬ 
tion  is  also  useful  in  the  problems  connected  with  selection  of  important 
discriminant  functions  in  the  area  of  classification.  Fisher  (1939), 

Hsu  (1939),  and  Roy  (1939)  have  independently  derived  the  joint  distri¬ 
bution  of  the  eigenvalues  of  the  MANOVA  matrix  in  the  central  case. 

Hsu  (1941)  derived  the  above  distribution  in  the  noncentral  case  when 
the  sample  size  tends  to  infinity  and  the  underlying  distribution  is 
multivariate  normal.  In  proving  the  above  result,  Hsu  assumed  that  the 
ratios  of  the  sample  sizes  of  the  groups  to  the  total  sample  size  tend 
to  constants  in  the  limiting  case.  In  this  paper,  we  extend, the  result 
of  Hsu  to  the  case  when  the  underlying  distribution  is  not  necessarily 
multivariate  normal. 

'  '  □ 


I 
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2.  PRELIMINARIES 

For  t  =  1,  2,...,k,  let  X'  =  (X, . .... X  )  be  distributed  with  mean 

~t  It  pt 

vector  =  (£.  )  and  covariance  matrix  Z  =  (a..).  Let  X  , — ,X 

t  1 1  p  t  tl  n  m  ^ 

be  a  sample  of  m^  independent  observations  on  X^.  Also,  let 

1  k 

=  N  N  =  V‘+mk 


■  <V  'I 

J  t=l 


(2.1) 


Unless  stated  otherwise,  we  assume  that  m^  =  q^m.  So,  N  =  qm  where 
q  =  q^+...+q^#  The  eigenvalues  of  H' Z  *  are  constants  independent  of  m. 
Let  the  rank  of  f  be  r  and  the  nonzero  eigenvalues  of  'fE  *  be  denoted  by 
\'.j>..._>V.  We  assume  that  these  eigenvalues  have  multiplicities,  i.e.. 


xi  “  xh  for  1  '  Vi  +  1*,**»V  h  =  lf  2 . v 


(2.2) 


where  aQ  *  0,  ah  *  a^  j  +  u^,  av  *  r*  Now»  let  4>j  >•  •  *2$ denote  the 

A  A  | 

eigenvalues  of  the  MANOVA  matrix  £  where  2^  =  min(p,  k  -  1)  , 

1  k  _ 

V  =  ±  l  m  (X  -  X..)(X  -  X..)’ 

t=l  ~ 


k  mt 


1  ■ »  J,  i,<?“ 


X  =  —  l  x„ .  , 
.t.  mt  .irt3 


X  )(X  -  X  )’  =  (a  ) 

~t.  ~t  2  -t.  lj 


1  k 

x. .  =  ^  y  >  xt.. 

N  t-i  j=rt3 


(2.3) 


Here  we  note  that  E(£)  =  —  —  k  I  and  E('f)  *  V  +  ^k  —  E. 

N  N 


In  the  sequel,  we  need  the  following  lemmas: 

Lemma  2.1  For  each  t,  t  =  1,  2,...,k,  let  :  p  *  1 ,  2  =  1 ,  2,..., 

be  a  sequence  of  i.i.d.  random  vectors  with  Ev^  =  0  an<^  with  I  as  the 

covariance  matrix.  Suppose  that  the  k  sequences  are  independent  and 

2 

that  for  each  t,  Efy^y^)  < 


00 


3 


If  m  =  q  m,  q  >  0,  t  =  1,  2,...,k  and  N  =  m.  +  - • •+m,  ,  then,  as  m  tends 

t  t  t  i  iC 

to  infinity,  the  k  vectors  (y  -  y.  .)  ,  t  =  1,  2,...,k  and  the  p  x  p 

L  ^  U  • 

matrix 


m 


^  \  l  XXtt  ■  -  z't ?  -  IP1 

VN  t=l  2=1  r 

converge  in  distribution  to  k  random  vectors  Z t  =  (z ^ ^ . .  .z  ) ' ,  t  *  1,  2,...,k 
and  a  p  x  p  symmetric  matrix  U  =  (u^ ^ )  ,  which  satisfy  the  following: 

1.  The  joint  distribution  of  ^ZjC»  j  =  1»  2 . p,  t  =  1,  2,...,k, 

u..,  1  f  i  j  j  £  pl  are  kp  +  — p  (p  +  l)-variate  normal , 
ij  2- 

2.  The  k  vectors  {Z^,...,Z^}  and  the  matrix  U  are  mutually 
independent , 

3.  EZ  =  0,  EU  =  0  and 


varZ 


t  =  (1  -  (qt/q))Ip,cov(Zt,  Zg)  =  -  (^q^  /q)Ip,  t  t  s 

cov(uii-  V  ■  {  J1it<E>,itiyjti 

cov<\vuW<  J.’tViV'ViV1 


h  *  h  or  J1  *  h>  1  -  V  V  h>  j2^  p' 


Here 


m 


,  m 
k  t 


t  2=1 


t=i  2=r 


i  i  ^  p  i  it 

Ft.  ■  -e  .1.  Xu-  ■  N  J,  Liti  ■  n  JjVt.- 


1^  denotes  the  p  y  p  identity  matrix  and  ^  is  the  j th  element  of  the 
p-vector  y  ,  and  q  =  q^+..*+qt. 

Proof :  The  proof  follows  by  application  of  central  limit  theorem  and 

by  direct  computations. 

I.emma  2.2  Let  X  ,  n  =  0 ,  1  ,  2 , . . . ,  be  a  sequence  of  random  p-vectors 
-  -n 

with  X^  -*■  Xq  in  distribution.  Then  we  can  find  a  probability  space 
('' ,  F,  P)  on  which  we  can  define  a  sequence  of  random  vectors  Xn> 
n  =  0,  1,  2,...,  such  that 


1.  X  and  X  are  identically  distributed 
~n  ~n 

2.  X^  -*•  Xq  pointwise. 

The  above  lemma  was  given  in  Skorokhod  (1956). 


L emma  2 . 3  Let  g^(x)  be  a  sequence  of  K-degree  polynomials  with  roots 
x|n\...,X^n^  for  each  n,  and  let  g(x)  be  a  k-degree  polynomial  with  roots 
Xj,...,X^,  k  _<  K.  If  gn(x)  g(x)  as  n  then  after  suitable 


rearrangement  of  x' 


(n)  ,  (n)  ,  „  .  , 

,x^  ,  we  have  x.  -*x.,j  =  l,  2,...,k  and 

K  3  J 


>,  j  =  k  +  1 ,. . . ,K. 


The  proof  of  the  above  lemma  is  given  in  Bai  (1984).  In  the  sequel. 


det  A  denotes  the  determinant  of  A. 
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3.  ASYMPTOTIC  JOINT  DISTRIBUTION  OF  THE  EIGENVALUES  OF 

THE  MANOVA  MATRIX 


In  this  section,  we  derive  the  asymptotic  joint  distribution  of  the 
eigenvalues  of  the  MANOVA  matrix  when  the  population  eigenvalues  have 


multiplicities.  Let 


n  =  E  £  and  n  =  E  ~ 


1  ^mt^t  "  h)  (nt  -  n)  '  =  E  ^  ||  ||  E 


which  has  \!,...,X*  as  its  nonzero  roots.  Consider  the  matrix 
1  r 


A  =  ( /q  1  /q  (n1  -  n)  ...  »/c[  / q  (n^  -  n)). 


It  is  evident  that  AA'  =  E  2  ||  il'.y  )|  E  ^  and  that  A(/q^/q  )  =  0. 

Let  B  ^  be  an  orthogonal  matrix  with  (  /q  ^  /  q  , . . . ,  /q^/ q  )’  as  its  last 
column,  and  A^  be  the  p  x  (k  -  1)  matrix  constructed  by  the  first  k  -  1 
columns  of  AB^,.  Noting  the  last  column  of  AB^  is  zero,  we  know  that 


AXA'  =  z'"2  ||  .  ||  Z' 


Xta =  z'2  -  5t> 


(3.1) 


1  t  l  kl 

It.  "  nT  J.Ztl.  and  '  N  J,Vt. 


t  ?,=  1 


k  mt 


=  t  J,  ‘  Zt^Vti  -  xij  ‘  V 


Hm  =  "  y..)..../!^  (yk  -  y..)) 


According  to  Lemma  2.1,  we  have 


(H  ,  Um> 

m  m 


in  dist. 


(H,  U) 


where  H  =  (Z  ,...,Z  ),  and  U  =  (u..)  are  mutually  independent,  and  are 

1  iC  1 J 

distributed  respectively  as  (Z  ,...,Z^)  and  U  defined  in  Lemma  2.1. 

-  g-1 

By  simple  computation,  we  find  that  the  eigenvalues  of  '-  .-  are 


the  same  as  the  solutions 

det  !|  ]  H  H* - —HA’  -  —AH'  +  AA’  -  $1  -  —  U  ||  =  0. 


(3.2) 


Let  G  be  the  p  x  (k  -  1)  matrix  constructed  bv  the  first  k  -  1 
m 

columns  of  H  B,.  Noting  that  the  last  column  of  H  B.  is  zero,  we  find  (3.2)  is 
ml  ml 

equivalent  to 


det  I)  ^  G  G’  -  —  G  A!  -  —  AG’  +  A.A!  -  M  -  —  U  j|  =  0  . 
N  m  m  m  1  1  m  11  ^  m 


(3.3) 


Note  that  the  p  x  (k  -  1)  elements  of  the  first  k  -  1  columns  of 
HBj  are  i.i.d.  N(0,l)’s  and  the  last  column  of  HB^  is  zero.  We  know 
that  the  elements  of  G^  converge  in  distribution  to  p  x  (k  -  1)  i.i.d, 

N (0 , 1 )  variables,  which  are  independent  of  U.  Let  be  the  p  x  (k  -  1) 
matrix  constructed  by  the  first  (k  -  1)  columns  of  HB^ . 

Recalling  that  AjA^  =  Z  *  [|  ^  ||  T.  ,  we  find  that  there  exists  a 

p  *  p  orthogonal  matrix  and  a  (k  -  1)  x  (k  -  1) 

Aj  =  B2  diagtVT^, . . .  0...01  B^ 

where  the  diag  K  >  J .  ,/r,  0...0]  denotes  a  p  *  (k  -  1)  matrix  whose 
first  r  diagonal  elements  are  Xj,...,X^  respectively,  and  the  remaining 
elements  are  all  zero. 


Let  W  =  BIG  BI  =  ||  w  m)||  ,  U  =  BlU  B„  =  ||  u.T'H  and  let  v  =  N"1  . 

IT!  2m3  1111  m  /  m  /  11  1  -i  1 1 


ij  11  m  2  m  2  "  ij 


Then  (3.3)  is  equivalent  to 


det  I!  v2W’W  -  vCfm)  +  D  +  vU  ||  =  0 
"mm  m  " 


(3.4) 


where 


D 


'  IT 

:  E1 


(m)  (m)  , 


c(m)  E(m) ' 
vv  v  | 

.  .  E(m)  0  1 

v 


.  (m)  _ 


V^h(wij’>  +  j  =  Vl  +  1 . *h  ’  h  U  2,-”,V’ 


X  w.  .  +  vi  w;"'  i  =  a,  ,  +  1 , .  . .  ,a,  ,  j  =  a  ,  +  1 , . . .  ,a 

g  ij  h  ji  n- 1  ’Ti  g-1  i 


1  <  h  4  g  <  v 


vT"  wfm) 

h  ij 


i  r  +  -1  *  j  - 


—  3.^.  i 


As  proved  earlier,  the  elements  of  H  and  of  U  tend  in  distribution 

mm 

to  that  of  H  and  U.  Hence  the  elements  of  W  and  of  U  will  tend  in  dis¬ 
til  m 

tribution  to  that  of  H  =  H  B’  *  ||  h^ [ {  and  of  U  *  B'  U  B?  =  |[  u^_.||  , 


2  2 


satisfying  that  H  and  U  are  independent  and  h^ '  s  are  i.i.d.  N(0,1) 
variables . 

According  to  Lemma  2.2,  without  loss  of  generality,  we  can  assume  that 
this  convergence  is  pointwise.  From  (3.4)  it  is  easily  seen  that  for  all 
ij  e  !2,  when  m  is  large  enough,  equation  (3.4)  has  p  roots,  and  they  are 


given  by  <j>,  >  <t<0>.  .  .><fc  . 

1  —  2—  —  p 

Write  U  =  II  V  ,  II  where  V  , 
11  gh"  gh 


lKuij)!l  ^  ^ •  >ag» 


j  =  +  l,...,ah,l_<g,  h_<v  +  l,  av+1  =  p  and  denote  by  Cgh>  Eg, 

the  limits  of  Egm^  ’  resPectively  and  take  the  variable  transformation 

<J>  *  Xj  +  ?v  in  (3.4).  Multiplying  by  —  the  first  p ^  rows  and  the  first 


pj  columns  of  the  determinant  on  the  left  hand  side  of  (3.4)  and  making 


m  tend  to  infinity,  we  obtain  the  following  in  the  limit: 


Let  the  p  roots  of  (3.4)  after  the  variable  transformation  be  denoted  bv 

.  fir 
’il 


’i™^  =  ^*i  ~  *  =  1*  2,...,p.  Then  we  know  that 


Ct.  i  -  1,2,. 


(m) 


i  =  3]  +  l,...,p,  where  ,  i  =  1,  2,...,a^  are  the  roots  of 

det  I!  -  C  +  Un  -  51  I!  =  0.  (3.5) 


1 


Similarly,  if  we  denote  by  5^  =  (<4.  -  A,  )  /-/v  the  roots  of  (3.4)  after 

In  1  h 


variable  transformation  £  -  (<{>  -  A  )  /V\F ,  we  can  prove  that 


+  °°  if  i  < 


Vi 


’ih 


h  if  1  ‘  Vi  +  '>• 

-  »  if  i  _>  +  1 ,  h 


.a- 


~  1,  2  , .  . .  ,v 


where  X^ ,  i  =  a^_^  +  l,...,a^  are  roots  of 


det  -  C,  ,  +  U,,  -  5 1 
1  hh  hh  u. 


=  0. 


(3.6) 


Finally,  if  we  take  the  variable  transformation  <J>  =  v  5  in  (3.4), 


1 


multiplying  by  —  the  last  p  -  r  rows  and  columns  of  the  determinantal  equation 


in  (3.4),  and  making  m  -*■  we  get 


det 


AT 

0 

e; 

1  U1 

i 

0 

...  A  I 

E' 

V  Uv 

V 

E, 

•  • .  E 

W-  51  ,, 

1 

V 

p-r 

=  0 


(3.7) 


where  W  is  the  (p  -  r)  *  (p  -  r)  right  lower  submatrix  of 


The  equation  (3.7)  is  equivalent  to 


9 


b 


, 


i 


det  ||  W  -  E  E!  -  "t“  E„E!  E  E'  -  a 

A1  l  1  a2  z  z 


If  we  write 


=  0. 


(3.8) 


B2H1B3  = 


w 

V 

w 

11 

•••  lk-1 

• 

• 

• 

• 

• 

• 

w 

w 

pi 

• • •  pk- 1 

where  the  w's  are  i.i.d.  N(0,1)  variables,  then  recalling  the  definition 
of  Ej,...,E  and  W,  we  find  that 


v;- 1' dijii  • 

1  V  J 


k-  1 


where  d..  =  w.,w.,,  i ,  j  =  r  4-  Hence  (3.8)  is  equivalent  to 

1J  V=r+1  U  ^ 


det(  l|  d.  .||  -  Cl  \  -  0. 

A'1  ij"  P~r/ 


(3.9) 


By  Lemma  2.3,  we  see  that 


r  (m)  _ 

t. .  .  .  =  N  cb . 

iv+1  x 


+  °°  if  i  <  r 


Ci  if  i  =  r  +  1 ,.  . .  ,p 


where  i  =  r  +  I,...,p  are  the  roots  of  (3.9).  Let  P.  ^  =  min(p,  k  -  I). 

Note  that  the  rank  of  ||  d_  ||  is  8.  j-r  with  probability  one.  Thus,  there 
are  surely  p  -  1.^  zero  roots  of  (3.9). 


Up  to  now,  we  have  proved  the  following  theorem. 

( 

1 


Theorem  4.1  When  m  tends  to  infinity,  then  c converge  in 


(  )  1 

distribution  to  ,  where  C  .  =  "^N  ($.  -  X ! )  ,  i  =  1,  2,...,r 

18-  l  li 

/  \  1 

and  t, ;  =  N<J» .  ,  i  =  r  +  1 , .  . . ,  8 ,  ,  and  for  n  =  1,  2,...,v,  c  +  1 , .  .  .  ,  C,  , 

111  Vi  A 

in  decreasing  order,  are  the  roots  of 


det  j  C  +  U.  -  Cl 
hh  hh  Vi 


=  0. 


(3.10) 


10 


S  are  Che  positive  roots  of 

r+ 1  /  j 

deC  li  11  diJ  11  ■  [Vr 


(3.11) 


C,  ,  —  j ;  ^  A,  ( w .  .  4-  w.  . )  I|  i,  j  !*•••*  3^  ^  1  *  2,..,,', 

hh  h  ii  li  ll  n-I  n 


“u  ■ ,  i  *• j  ■ r + 1 . »• 

J  i=r+l  J 


{w  ,  i  =  1,  2  , .  . .  ,p  ,  j  =  1,  2 , . . .  ,k  -  1}  are  i.i.d.  N(0,1)  variable, and 
are  independent  of  (U^,  h  =  1,  =  jj  u„  j|  >  i  >  j  =  +  1 

satisfying 


1)  [u..,  1  <  i  <  i  <  p  with  u  =  u..  are  —  p(p+l) -variate  normal 

ij  ~  -  ij  ji  2 

2)  Eu  =  0,  i ,  j  =  1  ,  2,.  . .  ,p 

k 

Eu. .u. .  =  -  l  q  (ER?  R2  -  1)  1  <  i,  j  <  p 
li  JJ  q  t“j  t  it  Jt  -  - 


Eu.  .  u.  .  =  -  Y  q  ER .  R.  R .  R,  t 

llx2  J1J2  q  t=l  C  V  V  V  V 


if  it  ^  i2  °r  ij  ^  i2- 

Here  R  are  the  elements  of 
it 


<h  -  “i> 


‘  ."k> 


R  .  •  •  •  R  . 
pi  pk 


Remark  3.1  The  admissibility  of  changing  -  in  (3.6)  into  in 

(3.10)  can  be  seen  from  the  symmetry  of  normality  of  the  entries  of 


C,  ,  and  the  independence  of  C,  ,  and  U. 
hh  K  hh 

Remark  3.2  In  Theorem  3.1,  it  is  easy  to  see  that  ( f; r+  j»  -  **»Cp  )  is 

independent  of  (r  , .  .  .  ,£  )  .  According  to  Hsu's  paper,  the  density  of  the 

distribution  of  ( (.  .,,...,5}  )  is  given  by 
r+ 1  **  j 


11 


-1 


D1  (Cr4-l . 2 


-  y  (p-r)  (k-l-r)  yUj-r))  ! 


i=r+l  j=i+l  1 


"1  /  "1 

n  (c.  -  On  n  c. 

i+l  3  )li=r+l 


>rr,jV4-|-i  +  l)| 

i-i  2  2  2  2'i 

5‘l  expr 1  J+.^j 


0  >  Oil- >  0, 
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where  =  max(p,  k  -  1),  =  niin(p,  k  -  1)  . 

It  is  evident  that  we  have  the  following  theorems. 
Theorem  3.2  Under  the  restrictions  of  Theorem  3.1  and  that 


ER.  R.  R.  R.  =0,  for  all  t  =  1,  2,...,k  and  at  least  one  among  i  i„j.j„ 

1 1  12t  121/ 

2  2 

is  not  equal  to  any  one  of  the  rest  and  that  ER^R.^  =  1,  for  all  t  =  1,  2,...,k 

and  i  ^  j,  the  v  +  1  random  vectors  (5,  , ... ),  (Ca  ......  )  .  . . (£ 

i  a.tl  a  ~  a  ..-rl 


1  a{  aj+l  a2  av-i+1  r 

(C^^...^  )  are  mutually  independent. 

Remark  3.3  If  for  each  t  =  1,  2,...,k,  has  an  isotropic  distribution 
with  mean  vector  y  and  covariance  matrix  E,  then  the  additional  condition 
of  Theorem  3.2  is  fulfilled. 
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Theorem  3.3  Under  the  restrictions  of  Theorem  3.2  and  that  ER.  =  3, 

-  it 


for  all  t  *  1,  2,...,k,  the  density  of  the  distribution  of 

t  y-  ..  \  I  /.  1  2  ,  .  .  ...  _  y  ' 


+1  >•••>£  )72Xh  +  4Xh  is  given  by  D(xa.  +1  ♦•••>*  ).  h  *  1 

h-1  n  Ti-1  h 


=  1 ,  2, ...  t  ,  where 


D(Xj . . .x^)  =  2 


L'.’w  {A  <2-22> 


(®  _>  2.  x2  •  '—xv  >  “  °°) 

When  the  underlying  distribution  is  multivariate  normal,  Hsu  (1940) 
proved  the  expression  (3.12)  and  (3.13).  In  proving  the  above  result,  Hsu  used 
Lemma  1  in  his  paper.  Recently,  one  of  the  authors  (Liang)  pointed  out, 
through  a  counterexample,  that  the  above  lemma  is  incorrect.  However, 

Bai  (1984)  had  very  recently  pointed  out  that  the  final  result  of  Hsu  is 


correct . 
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4.  APPLICATIONS 


In  a  number  of  situations,  it  is  of  interest  to  test  the  hypothesis 
Hj  :  =...=  Let  and  S w  respectively  denote  the  between  group 

and  within  group  sums  of  squares  and  cross  products  (SP)  matrices  respec- 

-  -  _1 
tivelv  where  SD  =  NT  and  SIt  =  NS.  The  nonzero  eigenvalues  of  S^S,.  are 
D  W  D  W 

given  by  $  _>.  .  •  It  is  known  that 

1  1 


(n  -  k  -  P  -  i)  -i  _  a'1 

(k  -  1)  E  SB  w  p  N(k  -  1) 


(4.1) 


The  problem  of  testing  for  the  rank  of  is  the  same  as  the  problem  of 

testing  for  the  number  of  significant  discriminant  functions.  A  dis¬ 
cussion  of  some  procedures  for  testing  for  the  rank  of  V,  was  given  in 
Krishnaiah  (1981).  We  will  discuss  as  to  how  the  results  of  Section  3 
of  this  paper  can  be  used  in  testing  for  the  rank  of  f.  when  the  under¬ 
lying  distribution  is  not  multivariate  normal  but  the  first  four  moments 
of  the  distribution  exist.  Unless  otherwise  stated,  we  assume  in  the 
sequel  that  (m^/m)  =  qt  for  t  =  1,  2,...,k  and  N  =  mq  tends  to  infinity. 

We  can  test  (e.g.,  see  Krishnaiah  (1981))  for  the  hypothesis 
by  using  certain  functions  q(  4^  ^ , . .  . ,  <^^ )  of  the  eigenvalues  of  TI 
where  denotes  the  hypothesis  that  the  rank  of  tl  *  is  r.  Fisher  (1938) 
suggested  using  T^  as  a  test  statistic  for  testing  where  T^  *  if r+j  +  *  •  .+4>? 
When  is  true,  the  joint  asymptotic  distribution  of  N ^  »•  •  •  is  the 

same  as  the  joint  distribution  of  the  eigenvalues  of  :  (e x  -  r)  x  (fc  -  r) 

where  S  is  distributed  as  the  central  Wishart  matrix  with  degrees  of 

freedom.  We  can  use  <fr+j  also  as  a  test  statistic.  The 

asymptotic  distribution  of  N4>r+j  ,  under  Hr>  is  the  same  as  the  distribution 
of  the  largest  eigenvalue  of  Sr>  This  distribution  has  been  tabulated 
(e.g.,  see  Krishnaiah  (1981)).  Suppose  the  null  hypothesis  is  not  true 


and  the  rank  of  is  r  +  s.  Then  ^+s+l~*  '  an°  *r+s  > 

Also,  if  are  distinct,  then 

r+i  r+s 


are  distributed  independently  as  normal  with  mean  zero  and  variance  one. 

So,  we  can  compute  the  asymptotic  power  function  of  the  test  based  upon 

<J>  , . .  If  the  eigenvalues  X'  .,..., X  ’  (t  <  s)  are  equal,  then  the  joint 

r+1  r-ri  r+t 

asymptotic  distribution  of 


is  the  same  as  the  joint  distribution  of  the  eigenvalues  of  the  Gaussian 
matrix  A  =  (a^)  :  (r  +  t)  x  (r  +  t)  where  the  elements  of  A  are  distributed 
independently  as  normal  with  means  zero,  varCa^)  =  2  and  var(a„)  =  1 
for  i  ^  j  .  The  distribution  of  the  largest  eigenvalue  of  this  matrix 
can  be  computed  by  using  the  method  discussed  in  Krishnaiah  and  Chang  (1971). 
Percentage  points  of  this  distribution  are  given  in  Krishnaiah  and 
Schuurmann  (1985).  So,  the  power  function  of  the  test  in  this  case  can 
be  computed. 

We  will  now  discuss  a  sequential  method  (see  Krishnaiah  (1981)) 
method  of  testing  for  the  rank  of  £2.  The  hypothesis  £2  =  0  is  accepted 


or  rejected  according  as 


where 


Pi*!  <  Col  |  £2  -  0]  -  (1  -  Oj). 
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If  n  =  0,  we  don’t  proceed  further.  If  C  =  0  is  rejected,  we  accept 


reject  according  as 


fc2  >  Ca2 


where 


P^2^Ca2  >  Cal;  V  =  °  -  V 


When  Hj  is  true,  <J>  and  <J> 2  are  distributed  independent  of  each  other 
asymptotically.  Also,  under  Hj ,  is  distributed  asymptotically  as 

the  largest  eigenvalue  of  :  (£^  -  1)  *  (£^  -  1)  where  is  distri¬ 
buted  as  the  central  Wishart  matrix  with  degrees  of  freedom  and 
ECSj)  =  2. ^ I •  So.  we  can  compute  a2  for  given  value  of  Ca2  and  vice 
versa.  If  is  accepted,  we  don’t  proceed  further.  Otherwise,  we 
accept  or  reject  according  as 


^3  >  Ca3 


where 


PI+3±Ca3 


h  -  Ca2 ’  H2] 


(1  -  a^)  . 


When  H2  is  true,  N<j>3  is  distributed  independent  of  <j>2  as  the  largest 
eigenvalue  of  S2  :  ^  -  2)  x  (j^  -  2)  where  S2  is  distributed  as  the 

central  Wishart  matrix  with  &2  degrees  of  freedom  and  E(S2)  =  £2I.  So, 
we  can  compute  C  ^  for  given  value  of  a ^  and  vice  versa.  This  procedure 
is  continued  until  a  decision  is  made  about  the  rank  of  fi. 
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